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Why the Effort? 


The trend in design indicates that future airplanes will be statically unstable in pitch, be more flexible 
than current aircraft, and require highly integrated, interdisciplinary, design methodologies [1]. Fighter 
aircraft will be more maneuverable and will use active flutter suppression. One application of active 
flutter suppression is to provide the required margin between the maximum attainable speed in a dive and 
the speed at which flutter occurs while also requiring open-loop flutter to be at or above the maximum 
dive speed. A more ambitious application of flutter suppression would be to allow the normal operating 
envelope to exceed the open-loop flutter speed. If active flutter suppression is to become part of the 
integrated flight control system, then an integrated modeling and simulation capability is required. This 
modeling and simulation capability would embrace traditional non-linear, rigid-body mechanics for 
aircraft and traditional linear aeroservoelastic dynamic models. In particular, a unified set of equations 
and notation should arise. 

A variety of programmatic responses arose from the concern that current modeling practices needed 
to be reexamined in light of anticipated applications to future aircraft. At the Langley Research Center a 
Functional Integration Technology (FIT) team was established to perform dynamics integration research 
using the F/A-18 as a focus vehicle. A central part of this effort has been the reexamination of the 
aeroelastic equations of motion for fixed-wing aircraft [2,3] and the development of a comprehensive 
simulation modeling capability [4]. At the Wright Research and Development Center, a 30-month 
contract was awarded to Lockheed to develop an aeroservoelastic analysis and design software package 
wherein the equations of motion are developed from first principles [5]. At the Air Force Office of 
Scientific Research a contract was let to Professor Luigi Morino to develop the equations of motion of a 
maneuvering, flexible airplane with minimal simplifying assumptions [6]. The Lockheed effort [5] adapts 
the method of hybrid coordinates used by Likens for space-craft applications [7] to the aircraft problem. 
Morino's approach [6] is very similar to the FIT effort [2,3] and does a nice job of incorporating the total 
vehicle rotational degrees-of-freedom in a Lagrangian framework by taking partial derivatives of kinetic 
and potential energy with respect to the entire direction cosine matrix. 


• OBSERVED TRENDS IN AIRPLANE DESIGN (FIGHTERS) 

- STATIC INSTABILITY IN PITCH 

- MORE FLEXIBLE 

* MORE MANEUVERABLE 

- HIGHLY INTEGRATED DESIGN 

- ACTIVE FLUTTER SUPPRESSION TO PROVIDE MARGIN 

. URGE TO "UNIFY" AND GENERALIZE NOTATION AND EQUATIONS 

- TRADITIONAL, RIGID-BODY AEROMECHANICS 

- TRADITIONAL, LINEAR ASE ANALYSIS 

• PROGRAMATIC RESPONSE - (1986-88) 

- LARC - FIT (Functional Integration Technology) 

- AFWAL - LOCKHEED ASE CONTRACT 

- USAF OFFICE OF SCIENTIFIC RESEARCH - L. MORINO 
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The Path Followed by FIT 


When the FIT team began its investigations, no references in the aircraft literature could be found 
wherein equations of motion for elastic airplanes were developed with what seemed to be sufficient 
generality. The usual developments made assumptions that resulted in the absence of any inertial 
coupling between angular rates and elastic deformation. Occasionally questions were raised as to the 
conditions under which these simplifying assumptions might not be justified [8]. The literature for 
spacecraft and rotorcraft was not examined initially. 

In the FIT team development, a Lagrangian approach was used to derive the equations of motion of 
an elastic airplane flexing about a rotating reference frame. As a result of retaining the coupling terms in 
the kinetic energy expression, non-linear terms representing inertial coupling between angular and elastic 
degrees-of-freedom were identified. Equations including these terms were implemented in a simulation 
model of an F/A-18 and a number of trajectories calculated to determine the effects of these coupling 
terms [2]. At the same time, a number of articles were appearing in the literature that examined the errors 
that can arise in predicting centrifugal stiffening when all the nonlinear terms are retained in the kinetic 
energy expression but not in the potential energy expression [9,10]. The result was that either zero 
stiffening or negative stiffening is predicted in cases where positive stiffening should result. This error 
was common to most of the multi-flexible body simulation codes then available and to the initial 
developments of the FIT team and Lockheed [2,5]. The approach adopted by the FIT team to correct the 
error in centrifugal effects was to augment the existing simulation structure with terms resulting from 
retaining non-linear strain/displacement terms in the potential energy expression and mapping their effect 
into the existing simulation structure. The FIT approach and its application to an example problem are 
described in detail in reference [3]. 


• PERFORMED LITERATURE SEARCH 

- VERY LITTLE IN AIRCRAFT LITERATURE ON MANEUVERING, FLEXIBLE 
STRUCTURES 

- SPACECRAFT AND ROTORCRAFT LITERATURE IGNORED 

• REDERIVED EQUATIONS OF MOTION 

- AIRPLANE FLEXING ABOUT A ROTATING FRAME 

- INERTIAL COUPLING BETWEEN ANGULAR MOTION AND ELASTIC MODES 
IDENTIFIED 

- APPLIED TO MODEL OF F18 TO DETERMINE EFFECTS - AIAA 87-2501-CP 

• THEORETICAL WEAKNESS ADDRESSED 

- CENTRIFUGAL EFFECTS - [Kane], [Eke] 

- APPROACH CHOSEN: NONLINEAR STRAIN/DISPLACEMENT - GEOMETRIC 
STIFFENING 

- APPLIED TO SIMPLE ROTATING STRUCTURE - AIAA 88-2232-CP 
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The Issues 


One is left with with two questions at this point: (1), Is geometrically-exact flexible body modeling 
theory required for the dynamic analysis of fixed-wing aircraft, even highly-augmented, unstable ones; 
and (2), how is such modeling best accomplished? Previous work [2] indicates that for most fixed-wing 
aircraft in a clean-wing configuration, geometrically exact theory is probably unnecessary. Exceptions 
may occur in cases of underwing-store and T-tail configurations. Nonlinear terms due to the complex 
geometry can assert themselves in the form of parametrically excited oscillations. The most compelling 
reason for geometrically-exact modeling theory may be just the urge to have a comprehensive theory that 
works for all cases. The difficulty is establishing a non-linear theory that is sufficiently correct for the 
airplane problem without resorting to full-blown computational structural dynamics. An example of the 
computational structural mechanics approach is the Large Angle Transient DYN amics (LATDYN) code 
[11] developed at Langley and the work of Belytschko and Hsieh on which LATDYN was based [12]. 
These methods [1 1,12] assign a separate reference frame to each finite element that translates and rotates 
("convects") with each finite-element This "convected" coordinate method has been applied to large 
deformation problems such as car collisions and antenna deployment and would accommodate the 
nonlinear rotational/elastic coupling of a typical airplane structure with ease. However, a theory based on 
assumed modes and correct to second order may be sufficient for airplane applications and would require 
less computer resources than an approach similar to that of [1 1] and [12]. 


(1) ARE GEOMETRICALLY-EXACT MODELING THEORIES REQUIRED 
FOR AIRPLANE DYNAMIC ANALYSIS? 

- SPACECRAFT AND HELICOPTERS - YES 

- MOST AIRPLANES - NO 

- REASONS FOR / POSSIBLE EXCEPTIONS: 

. INERTIAL FORCES APPROACH AERODYNAMIC FORCES 

• COMPLEX GEOMETRIES - STORES 

• THEORETICAL PURITY 

(2) IF SO, WHAT IS THE BEST WAY TO IMPLEMENT SUCH THEORY? 

- MUST PROPERLY MODEL CENTRIFUGAL EFFECTS OR LEAVE THEM OUT 

- MANY POSSIBLE APPROACHES 
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Outline 


The remainder of this paper is organized as follows. First the equations of motion are briefly 
described. Since an energy approach was taken in the development of the equations of motion [2,3], 
expressions for kinetic and potential energy are defined. The differences between the FIT model and the 
more typical aircraft aeroelastic equations [8,13] are explained. Prior to defining the the potential energy, 
a simple example [10] is presented to illustrate the notion of "geometric" stiffness. The higher order 
terms in the FIT potential energy expression are explained in light of the simple example. Once the 
equations are established, the effects of the including the nonlinear inertial coupling terms in the 
simulation model of an F/A-18 are presented. Time responses due to a high-authority roll command are 
compared for the following cases: (1), additional terms included; and (2), additional terms ignored. 
Finally, conclusions and recommendations are offered. 


• EQUATIONS OF MOTION 

- KINETIC ENERGY 

- EXAMPLE OF GEOMETRIC STIFFNESS 

- POTENTIAL ENERGY 

• EFFECT OF INERTIAL LOADING ■ F/A-18 

• CONCLUSIONS 
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Modeling Assumptions 


The operating assumptions in the equations of motion are listed below. Assumptions (2) through 
(4) are typical for airplane applications [8,13]. Assumptions (5) through (7) are atypical and lead to the 
differences between the FIT model and the more typical aeroelastic equations of motion [8,13]. 

The structural finite-element model obtained for the F/A-18 was a lumped mass model, which 
provided the principal motivation for assumption (1). For a continuum model, summations over the 
lumped masses are replaced with integrations over the entire airplane. The finite-element model had both 
lumped masses and lumped inertial quantities and both were utilized in the calculations. Assumption (2) 
and (3) are consistent with each other. Assumption (2) leads to a generalized Hooke's law and 
assumption (3) allows the superposition of deformation modes. Assumption (4) reflects the fact that 
gravity gradients are only of concern in spacecraft dynamics. 

Assumption (5) acknowledges that while deformation is assumed to be small, total vehicle angular 
rate may not be small. The result is that products of total angular rate and deformation rates are retained in 
the kinetic energy expression. The effect of assumption (6) is that the term resulting from summing 
(integrating) the cross product of deformation with deformation rates over the total vehicle is retained in 
the kinetic energy expression. Assumption (7) recognizes that a transverse deformation of a beam will 
result in axial strain. This effect becomes critical in correctly predicting centrifugal effects and is best 
explained in the simple example that appears later. Apart from assumptions (5), (6), and (7), the 
following development closely parallels other developments of the aeroelastic equations [8,13]. 


(1) LUMPED MASS APPROXIMATION 

(2) LINEAR STRESS/STRAIN 

(3) DEFORMATION IS SMALL - SUPERPOSITION 

(4) GRAVITY CONSTANT OVER THE AIRPLANE 

FOLLOWING ATYPICAL AEROELASTIC ASSUMPTIONS WERE MADE 


(5) PRODUCTS OF ROTATION RATE AND DEFORMATION 
ARE NOT NEGLIGIBLE 

(6) ELASTIC DISPLACEMENT AND ELASTIC VELOCITIES 
MAY NOT BE PARALLEL 

(7) NON LINEAR STRAIN / DISPLACEMENT 
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Symbols and Definitions 


For the purposes of calculating kinetic energy, the total vehicle is viewed as a collection of small 
rigid bodies centered at the finite-element node locations. The lumped mass (dm), and the lumped inertia 

th 

[dl]. associated with the i— node can be interpreted as the result of performing integrations over the 

volume of the i^ rigid body using s as a variable of integration and the mass density, a, as a weighting 

function. The i^ 1 rigid body undergoes a translational deformation dj and a rotational deformation (^. 

The assumption of small deformations allows the rotation to be described as a vector and implemented as 

a cross product. The vector r locates the undeformed position of the i— rigid body in the vehicle body 

frame. The origin of the body frame is at the center of mass of the total vehicle when the vehicle is in the 

undeformed configuration. The vector R locates an arbitrary point of the i^ 1 rigid body in the inertial 
frame. 
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Total Kinetic Energy 


The total kinetic energy is calculated in two steps. First the kinetic energy of the i— rigid body 
(lumped mass) is found by performing the volume integration indicated in the brackets. The inertial 

velocity of an interior point is squared, multiplied by the mass density, a, and integrated over the small 
rigid body. All volume integral expressions involving the variable of integration $ can be resolved into 
the "known" parameters, (dm), and [dl]., that are defined in the previous figure ..A summation is then 

performed over the small rigid bodies indexed by i. If the deformations are described as a sum of spatial 
functions /mode shapes) and time-dependent generalized coordinates, a separation of variables is achieved 
The kinetic energy becomes a summation of terms where each term is a product involving time varying 
coordinates and constant mass or length terms resulting from the summation (integration) over the total 
vehicle. 


. INTEGRATE OVER THE IDEALIZED LUMPED MASS TO GET 
KINETIC ENERGY FOR EACH LUMPED MASS 

• SUM THE KINETIC ENERGIES ASSOCIATED WITH EACH 
LUMPED MASS 

t = fz { /// o(-£a)-(Jfi)<is } 

» OVER I UI I 


WHERE 

R, — P + r + d + s + 0 * § 


i = 

a. = sm 1 

i 1 


dt T 


TRANSLATIONAL DEFORMATION 
IN MODAL COMPONENTS 

ROTATIONAL DEFORMATION 
IN MODAL COMPONENTS 

TIME RATE OF CHANGE WITH 
RESPECT TO THE INERTIAL FRAME 
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Kinetic Energy in Modal Components 


The first three terms of the expression defining kinetic energy, shown below , are found 
in the standard aeroelastic equations of motion developments [8,13] and represent "rigid" translational 
energy, "rigid" rotational energy, and elastic kinetic energy, respectively. The term [J q ] is the inertia 

dyadic of the total undeformed airplane expressed in body-frame components and has units of mass- 
2 

length . The fourth line describes coupling between translational, rotational, and elastic momenta. If the 
assumed mode shapes are the modes of free vibration of an unrestrained structure, they satisfy the first- 
order mean axis conditions [8,13] and the terms, &• and h- are zero. The term a. is simply the location of 

th J J J 

the center of mass of the mode shape in the body frame and has units of length. The term h. is the first 

th ^ 

moment of the j— mode shape in the body frame using mass as the weighting function and has units of 
2 

mass-length . The terms in the dashed-line box result from assumptions (5) and (6). The terms [J]j and 
[JJjk are the first and second partial derivatives of the inertia dyadic matrix with respect to elastic modes j 
and k. The term results from summing (integrating) the cross product of mode shape j with mode k 

over the vehicle and has units of mass-length^. For a more complete explanation, see reference 2. 


T H y rnV-V • RIGID BODY TRANSLATIONAL KINETIC ENERGY 


+ y CO*[J 0 ] * RIGID BODY ROTATIONAL KINETIC ENERGY 

• « 

+ yM T| J 11 • MODAL ELASTIC KINETIC ENERGY 

jk 


mV»a.jV + + mV.coxajV 


• ZERO FOR UNRESTRAINED 
MODES 



NONLINEARITIES 
RESULTING FROM 
ASSUMPTIONS (5), (6) 
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Illustration of Geometric Stiffness 

The example shown below is taken from reference [10] and provides a simple paradigm of the beam 
problem. Pictured is a simple 3 degree-of-freedom, planar system. A mass, m, is located at the 
outermost point and is the only mass in the system. The structure rotates freely about the point P. The 

angle between rod a and an inertial reference is given by \| /. There is a torsional spring, k 0 , and a linear 

spring, k d - The deflection of the torsional spring is given by 0 and the linear spring by d = r5, where 5 is 

a nondimensional deflection. Another coordinate system (x,y) is given by an axis system located at the 
zero-strain position of the mass, m. The (x,y) coordinates are analogous to those typically used in the 

beam problem. The coordinates (^,rj) are the non-dimensional forms of (x,y). Both the (\|/,5,9) and the 

(V&Tl) coordinate systems are equally valid for describing this system. The (\jr,5,0) system leads to a 

particularly simple expression for strain energy. The reason is that a change in 5 produces only a linear 

distortion in the spring, kg, and similarly for 0. Thus the strain energy, U, is given by, 

U = (1/2) { k 0 0 2 + k g 8 2 } 

where k g = k d r . The (\|/,£,T|) system produces a complex expression for strain energy. A change in r| 

produces nonlinear changes in both springs and similarly for t, if there is some deflection in r\. Using the 
relations, 

6 = -1 + {T1 2 + (1+£) 2 } 1/2 and 0 = Tan'^H^)}, 
one gets for the strain energy, U, in terms of % and rj, 

U = (1/2) { kg{ -1 + [T | 2 + (1+^) 2 ] 1/2 } 2 + k 0 { Tan" 1 [T|/(l+^)] } 2 ). 



. TWO CHOICES FOR COORDINATES: (y , 8 , 0 ) OR ( ^ , r) ) 

• ( 8 , 0 ) -> SIMPLE EXPRESSION FOR STRAIN ENERGY 

• (^,q ) -> COMPLEX STRAIN ENERGY EXPRESSION 

• PARADIGM OF BEAM PROBLEM 
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Comparison of Resulting Linear Models 

If the full nonlinear equations are derived using the (\j/,S,0) coordinates, the \|/ equation removed, 

and the quantities y and \| r are treated as parameters, the remaining two equations can be linearized in 0 
and 8 and the resulting correct linear equations are given in case 1 below. Suppose the procedure is 
repeated for the (\|f,^,r|) coordinates except that a linear approximation for strain energy given by, 

U ~ (1/2) { k 0 ii 2 + k^ 2 }, 

is used, then the resulting equations are given as case 2 below [10]. The only differences occur in the 

stiffness matrix. The foremost difference is that a de-stiffening result is predicted in case 2 for the r| 
degree-of-fieedom due to the spin rate when a stiffening effect should be predicted. 
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Potential Energy 


Returning to the equations of motion development for an airplane, the total potential energy is 
composed of a component due to the work performed by strain and a component due to the work 
performed against gravity. The gravity component is straightforward since gravity is assumed constant 
over the airplane. If a second order expression is used for describing strain as a function of displacement 
and a linear stress / strain relationship exists, then a fourth-order strain energy expression results [3]. The 
strain energy expression, however, is accurate only to third order, so the fourth-order term is dropped. A 
third-order term in the strain energy expression, shown in modal components below, would lead to 
second-order stiffness terms in the final equations. If the third-order strain-energy term is left in this 
form, then modes with significant axial displacement need to be included in the dynamic model. As seen 
from the simple example, coupling occurs between axial and transverse displacement in the strain energy 
expression when beam-like coordinates are used. Since, for airplanes, axial modes are typically much 
higher frequency than transverse modes, the effects of the axial modes can be residualized. This 
residualization is accomplished in the FIT framework by solving for the elastic displacements that result 
from combinations of unit values of rotational velocity about the airplane body-frame axes. Thus for a 
unit roll rate, p, the steady state deflections are calculated for the full finite-element model. These "static" 
deflections are combined with the third order stiffness tensor to produce an increment to the basic 

stiffness matrix that is appropriate for unit roll rate. Thus the j,k^ element of an incremental stiffness 
matrix, [AK], is given by, 


AK jk = AK jkp<s • 

The j,k*k element of [AK] becomes one entry into the 3x3 matrix [J]^. Because of symmetry, the 

matrix, [J]^, has 6 free parameters, and each is calculated by repeating this process for different 

combinations of unit roll rate, p, pitch rate, q, and yaw rate, r. This process is discussed in greater detail 
in reference [3]. 


POTENTIAL ENERGY h STRAIN + GRAVITY CONTRIBUTION 

U = U + U 

s g 

STRAIN ENERGY IN TERMS OF MODAL STIFFNESS MATRICES 

U s = (1/2) [K + aK% p + aV r, P if ] r, J n k 

jk jkp jkpq 

= (i/2) [ - w-{l[J]^}*ffi]ri j ri k 

HIGHER ORDER TERMS RESULT FROM NONLINEAR STRAIN / 
DISPLACEMENT EQUATIONS 
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Final Equations - Unrestrained Modes 

The final equations in modal form are shown below. These equations apply to the case where the 
assumed modes satisfy the first-order mean axis conditions, otherwise additional coupling is present. 
The terms that are inside the dashed boxes result from assumptions (5), (6), and (7). The terms outside 
the boxes are equivalent to the equations seen in the more traditional approaches [8,13]. 


TRANSLATIONAL MOMENTUM 


m V + m coxV = F + mg 


ROTATIONAL MOMENTUM 




+ t> V i] \ + ax [J]m 

jk 



= L 


ELASTIC MODE j 


M_ if 
jk 


- m-ii n k + Kn k > aa-u rj k - j- a{ [J] . + [J] n + [J]% }a 

jk jk jk * J jk jk 


= Q. 


TOTAL INERTIA M ATRI X IN THE BODY FRAME 
[J] = [Jj|+ [J] j r| J + y { [J] + [J]p ri J r| k 
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Inertial Effects - Time History Responses 

In assessing the effects of the additional terms in the FIT equations of motion on a predicted 
response, a variety of cases were calculated. While all three axes were examined, the most interesting 
responses involved roll rate. The responses shown here were presented at the 1987 Flight Simulation 
Technologies Conference [2] and were generated prior to the incorporation of the geometric stiffness 

terms, [J]^, in the simulation model. Preliminary runs made subsequently, but not presented in this 

paper, suggest that the effect of the additional stiffness terms is small for the angular rates considered and 
that qualitative conclusions drawn from the data presented in reference [2] are still valid. 

The time responses were generated by injecting a roll command doublet at the actuator input. A 
combination of aileron and stabilator was used. The initial conditions were straight and level flight at 
Mach .7 at sea level. Responses were generated with ("terms on") and without ("terms off) the 
additional angular/elastic coupling terms that are part of the FIT equations of motion. 


. GOAL: ACHIEVE SUFFICIENT ANGULAR RATE TO EXCITE 
INERTIAL RESPONSE 

• ROLL DOUBLET - 1 SECOND EACH SIDE 

• 20 DEG AILERON / 10 DEG DIFFERENTIAL STABILATOR 

• INITIAL CONDITIONS 

- MACH = .7 

- SEA LEVEL 

- STRAIGHT AND LEVEL 1 G TRIM 

• COMPARE RESPONSE WITH AND WITHOUT ANGULAR/ELASTIC 
INERTIAL COUPLING 

- "TERMS ON" • FIT MODEL 

- ’TERMS OFF" - TYPICAL ASE MODEL 


828 



Lateral / Anti-Symmetric Response 

For the anti-symmetric responses to the roll command, there was no discemable difference between 
"terms on" and "terms off' responses. The responses shown below, roll rate and the third anti-symmetric 
mode, are typical. The third anti-symmetric mode is characterized as a wing first-torsion mode with 
significant missile pitch. The finite-element model had tip missiles. The reason that no difference 
occurred is that for anti-symmetric modes, the elastic modes are excited by the terms pq, pr, and qr in the 
angular/elastic coupling terms. Since the roll maneuver remained well-coordinated, the pitch rate, q, and 
the yaw rate, r, remained small. Thus the coupling terms, pq, pr, and qr, remained small. 


ROLL RATE 
RAD/SEC 



TERMS ON 

TERMS OFF 


MISSILE PITCH/ 
WING 

FIRST-TORSION 
ASM 3 
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Symmetric Response - Set 1 


The inertial coupling terms made no effect on rigid body symmetric responses. The angle-of-attack 
response is shown below. The angle-of-attack at time equal zero is due to the 1-g trim. The second 
response shown below is the first symmetric mode, wing first-bending. The y-axis indicates the 
deflection of the mode in feet measured at the point of maximum deflection, presumably at the wing tip. 
The mode is positive for tip-down deflection, so the mode is displaced about .25 feet tip-up at time zero. 
While the angle-of-attack is clearly the principal driver of symmetric wing first-bending in this maneuver, 
a discemable difference has occurred between the "terms on" and the "terms off responses. Only two 
symmetric modes showed more difference in "terms on/off' responses than the first symmetric mode and 
these two modes are shown in the figure on the next page. 


5 

° 0 

DEG 

-5 



TERMS ON 
TERMS OFF 


WING 
FIRST- 
BENDING 
SYM 1 



TIME SECONDS 
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Symmetric Response - Set 2 


Symmetric modes 4 and 7 were the only modes to show significant inertially induced response. 
Symmetric mode 4 is characterized by tip-missile yawing together with fin-bending and symmetric mode 
7 is a wing first-bending in the fore/aft plane. The fact that roll rate squared is the principal driver is seen 
clearly in the figure. The y-axis for the elastic mode responses shown below represents strain energy 
absorbed normalized by that of the first symmetric mode. In other words, a unit deflection in modes 4 or 
7 represents the same strain energy absorbed as would occur with a deflection of one foot in the first 
symmetric mode. Since the wing is swept backwards, a forward deflection in the the wing results in 
outboard movement of mass and is therefore excited by roll rate. A natural question to ask is if these sort 
of responses can be predicted from merely examining the parameters of the simulation model. 


SQUARE OF 
ROLL RATE 


2 2 
RAD / SEC 


MISSILE YAW / 
FIN BENDING 
SYM 4 


WING IN-PLANE 
LATERAL 
FIRST-BENDING 
(FORE/AFT) 

SYM 7 




0.1 

0 

- 0.1 



TERMS ON 
TERMS OFF 


0 


1 

TIME 


2 

SECONDS 
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Modal Sensitivity Parameter 


As part of determining when inertial coupling might be important in an analysis, the following 
simple parameter can be calculated. The parameter Rj(fij) defined below is simply a first order 
approximation to the steady- state response of elastic mode j to a constant angular rate represented by the 

angular velocity vector, <&. The y constant is used to scale the responses so that for different modes j, 
identical Rj's represent the same strain energy. 




f { ffl T [J]. a } 

* J 


K 

j 


• MODAL DISPLACEMENT NORMALIZED BY STRAIN ENERGY 

• APPROXIMATES RESPONSE OF MODE "J" TO CENTRIFUGAL 
LOADING FOR A GIVEN ROTATIONAL VELOCITY 

. MAX DESIGN ROTATIONAL RATES ARE LIKELY INPUTS 
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Modal Sensitivity at Max Roll 


The parameter Rj(ffi) is shown below for each of the 20 elastic modes included in the simulation 
model. The input angular rates correspond to the max roll rate achieved at about .9 second into the 
maneuver. The angular velocity vector © = [p,q,r] in body-frame components. The fourth and seventh 
symmetric modes are clearly singled out by the Rj(ffl) parameter. Again, the units of Rj(fQ) are strain 
energy normalized to the first symmetric elastic mode. 


R.(<a) 

3 ' 7 

MODE SYMMETRIC ANTISYMMETRIC 


1 -.0283 .0008 

2 .0048 -.0120 

3 .0064 .0004 

4 | .0739 | .0039 

5 .0075 -.0001 

6 -.0028 -.0090 

7 [ 0899 1 -.0016 

8 -.0115 -.0006 

9 -.0043 .0018 

10 -.0027 -.0012 


W = [ 4.7 , .17 , .03 ] (RAD/SEC) 
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Illustration of R Parameter 


The amplitudes of the response of symmetric modes 4 and 7 arc essentially equal to the 
values of the R parameter calculated in the previous figure. This parameter, which is a linear 
approximation to steady-state response, is far from the end of the story. In the case of the 
FIT simulation model, even though symmetric modes 4 and 7 were excited by inertial effects, 
these modes are essentially decoupled from the rest of the dynamic model. This occurred 
because both modes 4 and 7 are dominated by in-plane bending of the wing lifting surface. A 
doublet-lattice code was used to calculate the generalized aerodynamic forces and in-plane 
motions produce no change in the normal washes induced at the 3/4 chord points of the 
aerodynamic boxes. Thus none of the other modes are significantly affected by symmetric 
modes 4 and 7. 

One can imagine other cases where inertially affected modes are more coupled to the rest 
of the system dynamics. One case is if such a mode contributes to a feedback signal. 

Another case might occur in an underwing store configuration. As the underwing stores were 
slung outboard by centrifugal forces, they would induce out-of-plane bending in the wings, 
the primary lifting surfaces. 
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Conclusion 


An integrated, nonlinear simulation model suitable for aeroelastic modeling of fixed-wing aircraft 
has been developed. While the author realizes that the subject of modeling rotating, elastic structures is 
not closed, it is believed that the equations of motion developed and applied herein are correct to second 
order and are suitable for use with typical aircraft structures. The equations are not suitable for large 
elastic deformation. In addition, the modeling framework generalizes both the methods and terminology 
of non-linear rigid-body airplane simulation and traditional linear aeroelastic modeling. 

Concerning the importance of angular/elastic inertial coupling in the dynamic analysis of fixed-wing 
aircraft, the following may be said. The rigorous inclusion of said coupling is not without peril and must 
be approached with care. In keeping with die same engineering judgment that guided the development of 
the traditional aeroelastic equations, the effect of non-linear inertial effects for most airplane applications is 
expected to be small. A parameter has been presented to help in the determination of when such effects 
are significant. The parameter does not tell die whole story, however, and modes flagged by the 
parameter as significant also need to be checked to see if the coupling is not a one-way path, i.e. the 
inertially affected modes can influence other modes. Classically, configurations where nonlinear inertial 
effects can come into play are characterized by complex geometries such as stores mounted under the 
wings or the presence of a T-tail. 


• INTEGRATED NONLINEAR MODEL DEVELOPED 

- CORRECT TO SECOND ORDER 

- SUITABLE FOR AIRPLANE STRUCTURES 

• GENERALIZES CONVENTIONAL ASE MODELS AND NONLINEAR 
RIGID-BODY MODELS 

• ANGULAR / ELASTIC INERTIAL COUPLING 

- RIGOROUS INCLUSION PROBLEMATIC 

- EFFECT NORMALLY SMALL FOR AIRPLANES 

- EXCEPTIONS CHARACTERIZED BY 

• R (CD) PARAMETER IS SIGNIFICANT FOR SOME MODE 
AND 

. AFFECTED MODE IS COUPLED TO THE REST OF THE MODEL 
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